Bistable defect structures in blue phase devices 
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Blue phases (BPs) are liquid crystals made up by networks of defects, or disclination lines. While 
existing phase diagrams show a striking variety of competing metastable topologies for these net- 
works, very little is known as to how to kinetically reach a target structure, or how to switch from 
one to the other, which is of paramount importance for devices. We theoretically identify two con- 
fined blue phase I systems in which by applying an appropriate series of electric field it is possible 
to select one of two bistable defect patterns. Our results may be used to realise new generation and 
fast switching energy-saving bistable devices in ultrathin surface treated BPI wafers. 
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Liquid crystals are materials with spontaneously bro- 
ken symmetry: as such they support topological defects 
of various kinds and may be used as testing grounds for 
theories on many areas of physics, e.g. the dynamics of 
cosmic strings and of vortex lines in superfluids and su- 
perconductors pQ. Blue phases (BPs) are a spectacular 
example of soft matter in which a network of such defects 
self assembles as a fully periodic cubic three-dimensional 
structure whose local director has a double-twist config- 
uration [2]. The defects are arranged in a regular fash- 
ion, with periodicity at length scales comparable with the 
wavelenght of light. This is the reason why BPs exhibit 
selective Bragg reflection in the range of visible light and 
can have interesting applications in fast light modula- 
tors [3] and tunable photonic crystals |T. BPs were ini- 
tially observed to be stable only in a narrow range of tem- 
perature ~ IK, but recently new compounds, stable in 
wider interval over ~ 60K (including room temperature), 
have been developed [5J [5] . These studies can be consid- 
ered with good reason as seeds for the recent fabrication 
of the first blue phase display device with fast switching 
times [7] . On the other hand BPs have attracted broader 
interest among physicists because double-twist cylinders 
represent a striking example of the so-called skyrmions, 
topological excitations encountered in nuclear physics [5], 
in spinor Bose- Einstein condensates [3] , and in ferromag- 
netic materials, where they have been observed [TO] . 

Theory and experiments have mapped out a number 
of phase diagrams reporting a remarkably varied range 
of competing metastable disclination networks [2J [TTj . 
However, there are only very few theoretical studies to 
date regarding the dynamics of the formation of BPs, 
and virtually none addresses how to switch from one 
state to another reversibly in a controlled way. Here we 
provide the first theoretical study to achieve this goal, 
and give two specific examples of a bistable BPI sys- 
tem, where an electric field along the appropriate di- 
rection can select one of two metastable arrangement of 
defects. Besides being of interest per se, as an intrigu- 



ing example of controlled defect dynamics, our discovery 
suggests specific novel geometries which may be used to 
realise bistable BP based devices. Devices like the ones 
we have found are important technologically as each of 
the bistable states is metastable and is retained after the 
field is switched off: this fact eliminates the need of hav- 
ing a constant electric field to keep the system in the 
"on" state, hence sharply reducing energy consumptions. 
Some examples of such systems are the zenithal bistable 
device (ZBD) [T2] and surface-stabilised cholesteric tex- 
ture (SSCT) devices [13]. In the ZBD case bistability 
relies on flexoelectricity and requires a periodic grating 
of the surfaces, whereas in the SSTC case it exploits the 
existence of focal conic defects, which compete against 
the defect-free helical texture when strong anchoring is 
used. Several other nematic multistable devices also ex- 
ist (see e.g. [H]). Our case is different as the switching 
occurs between two disclination networks. BP devices 
such as the one we suggest should also share the faster 
switching times made possible by BP cells. At the same 
time, simulations are key to clarify the dynamics of the 
defect reorganisation underpinning bistability [T2j [14] . 

In this work we propose two novel BPI (blue phase 
with 0% symmetry) geometries, for which it is possible 
to design a simple switching on-off schedule that leads 
to a bistable defect network. In both cases, bistability 
is present for a sizeable range of parameters (quantified 
below) and is characterized by the presence of two states, 
which are both metastable when the field is off and have 
a fully 3D defect structure. The first device we pro- 
pose relies on surface memory effects which retain the 
BPI structure on the boundaries, and the two compet- 
ing states are the standard bulk BPI network and an- 
other arrangement, which is only slightly higher in free 
energy. The second device exploits the additional frus- 
tration driven by strong homeotropic anchoring in a con- 
fined thin sample, which conflicts with the bulk ordering. 
We find that this leads to the appearance of multiple de- 
fect states, similar to the ones observed in P~5](a), which, 
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importantly, can be uniquely selected by the application 
of an electric field along an appropriate direction. Our 
devices are entirely different from the BP bistable de- 
vice described in [IB]. Firstly, in that device the dielec- 
tric anisotropy is negative, while in our case it is posi- 
tive. Secondly, the switching in [TBJ is between two field- 
induced metastable states, which are selected by different 
values of an applied voltage, whereas in our case the elec- 
tric fields is switched on and off along different directions. 

The physics of a BP device may be described by 
a Landau-de Gennes free energy written in terms of 
the tensor order parameter Q a p [2]. This comprises 
a bulk term f b = f(l - %)Q 2 a0 - ^Q a0 Q M Q ia + 
^fiQlp) 2 , and a distortion term f d = f{d p Q af3 ) 2 + 
^■(caCsd^Qsp + ZqoQap) 2 , where K is the elastic constant 
and the pitch of the cholesteric liquid crystal is given by 
p = 27t/q . Aq is a constant, 7 controls the magnitude 
of the ordering and e a( ^s is the Levi-Civita antisymmetric 
third-rank tensor. The constants Aq, K and 7 are cho- 
sen in order to be in the appropriate region of the phase 
diagram [II]- We have set A a = 0.0034, K = 0.005, 
7 = 3.775 for which BPI is stable. This choice may be 
mapped to a rotational viscosity equal to 1 Poise, (Frank) 
elastic constants equal to ~ 30 pN and a blue phase 
periodicity equal to 300 nm. Moreover one space and 
time unit correspond to 0.0125 /zm and 0.0013 fis. The 
contribution of the boundaries is encoded in a surface 
free-energy term, given by f s = l/2W (Q a p ~ Q%) 2 , 
where Wo (larger than 0.01 N/m in our simulations) 
controls the anchoring strength and is the order 

parameter at the surfaces. The interaction with an 
external electric field E is accounted for by an addi- 
EaQapEp, where e a (here taken >0) 
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tional term, 

is the dielectric anisotropy of the material. The elec- 
tric field may be quantified via the dimensionless num- 
ber e 2 = 27e a E 2 c /(32TrA j). By assuming a dielectric 
anisotropy of ~ 100 [17], one finds that an electric field 
of 1 — 10 V/fim corresponds to e 2 ~ 0.003 — 0.3. The sim- 
ulation domains along the three coordinate axes are L x , 
L y and L z respectively, with periodic boundaries along 
x and y. All our results are performed on lattices of di- 
mensions 32 x 32 x 32. The equation of motion for Q is 
D t Q = TH, where T is a collective rotational diffusion 
constant and D t is the material derivative for rod-like 
molecules 18 . The molecular field H is the thermo- 
dynamic force that drives the system towards the equi- 
librium [18]. The velocity field obeys a Navier-Stokes 
equation with a stress tensor generalised to describe LC 
hydrodynamics [18] . The interplay between the velocity 
field and the order parameter is referred to as backflow. 
To solve the equations of motion, we used a hybrid lattice 
Boltzmann algorithm, as in |14j . 

The usual zero-field defect structure in the BPI cell is 
shown in Fig. la. We first focus on the case in which 
the director field is fixed at the top (z = L, being L the 





FIG. 1: Equilibrium disclination networks for BPI when a) 
the structure is strongly anchored at both walls, and when b) 
there is homeotropic anchoring at both walls. 



sample size) and bottom (z = 0) boundaries of the cu- 
bic cell to its stable structure in absence of a field. This 
may be achieved by pinning due to impurities or through 
surface memory effects [19]. Fig. 2a-d shows the evolu- 
tion of the defect network in response to an electric field 
switched on and off first along x, then along y. During 
the cycle the defect dynamics is highly non-trivial and is 
characterized by a complex reorganization driven by the 
electric field, and depending on its direction. In particu- 
lar, after the application of the field along x, the defects 
twist and bend and eventually recombine, annihilating 
in the bulk and forming defect arcs pinned at the walls 
(ti — ti). The voltage (AV y = E y L y ) is chosen strong 
enough to break the defect structure in the bulk. When 
the field is switched off, some disclinations migrate back 
to the centre of the cell, giving rise to connected bound- 
ary arcs that span the entire device (t 5 — 1 8 ). The system 
gets stuck into a stationary metastable state (tg), whose 
free energy is different from the one of the zero-field state 
(see inset of Fig. [2j top panel) . Due to the anisotropy of 
the BPs, it is now interesting to ask what is the defect 
dynamics if the electric field acts on a different direc- 
tion. Starting from the configuration at time tg we have 
switched on an electric field along the y direction, with 
voltage equal to the one chosen for the field in x direc- 
tion. The defects now start to deform, bending in the 
direction of the field (£9), and then move towards the 
bulk where they met and form unstable branching points 
that annihilate. This second field-induced state (£12) is 
again (as the one at time £4) characterized by boundary 
arcs. By switching off the field, the system reorganizes as 
the boundaries force the disclinations towards the bulk: 
surprisingly, the dynamics now leads back to the starting 
bulk BPI (Fig. la). Bistability is affected by the voltages 
used in the switching. We have simulated the same de- 
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vice for e 2 ~ 0.2, 0.3 and found that it is bistabile in both 
cases. On the other hand, bistability is lost for e 2 < 0.2 
and the switching dynamics is now characterized by a 
reversible cycle in which a zero-field state (different from 
the bulk BPI structure, Fig la), obtained after the ap- 
plication of the electric field along x, is recovered even 
after switching along y. Note also that the results in Fig. 
2 refer to BP devices with e a > 0: intriguingly, in the 
case of negative dielectric constant the device can never 
get back to the bulk BPI state. This may be due to the 
fact that in the presence of a field the director can form 
a helix along the field, and this is not different enough in 
free energy from the BP state to drive a reorganisation 
of the unit cell once the field is removed. To check the 
robustness of our results, we have also simulated the de- 
vice with different values of A and 7 (e.g. A = 0.0012 
and 7 = 5.085), again leading to a stable BPI - even 
in this different region of phase space we find that the 
device is bistable. Finally, in Fig. 2, top panel, we also 
compare the free energies as a function of time both in 
the presence and absence of hydrodynamics. The effect 
of backflow is to slightly speed up all relaxations, without 
deep consequences on the dynamic evolution. 

We now turn to the case of different boundary condi- 
tions. Fixed boundary conditions (FBCs), as the ones we 
have just employed, have often been used in cholesterics 
and BPs [15] . However in practice it is more feasible to 
control the alignment of the director field at the bound- 
aries by chemical treatment or rubbing. The former tech- 
nique typically is used to enforce homeotropic (or perpen- 
dicular) anchoring of the director at the wall, the latter 
yields homogeneous (or parallel) anchoring. Therefore it 
is interesting to know how these boundaries affect the 
switching dynamics of a BP device, and, importantly, 
whether they can still lead to a bistable behaviour. In 
what follows we focus on a BPI device with homeotropic 
anchoring at both walls. To the best of our knowledge 
there is only one study about the equilibrium disclina- 
tion configuration in confined BPI cell with homeotropic 
anchoring [15] (a) while nothing is known about electric 
field effects in these thin cells. As we shall see, by using 
the same electric field dynamical schedule adopted for the 
cell with FBCs, it is still possible to realise a switchable 
bistable device, albeit with some key conceptual differ- 
ences with respect to the device in Fig. la. The first one 
is that homeotropic anchoring affects the zero-field equi- 
librium defect structure (Fig.[T]b)), now characterized by 
bent disclinations depinned from the walls - the original 
topology of BPI being preserved. The second important 
difference is that the change of the boundary conditions 
affects the dynamic evolution under an electric field: the 
minimum value required to destabilise the initial topol- 
ogy is now e 2 ~ 0.05. This decrease in the threshold is 
because the wall now pins the zero-field structure more 
weakly than in the case of Fig. 1. In spite of these differ- 
ences, there is a range in which the device shows bistabil- 




FIG. 2: Free energy evolution (top panel) of the BPI device 
with FBCs, with (red) and without (black) backflow. Param- 
eters were A = 0.0034, K = 0.005, 7 = 3.775, e 2 ~ 0.2 
(corresponding to a field of ~ 8 V/fim). The evolution of the 
defects under an applied field is reported in the bottom panel. 
The step function in the top panel is 7^ when the field is 
on. From t = up to t = 35000 the system relaxes to the 
stable BPI configuration in Fig. la. An electric field along 
x is switched on at t = 35 x 10 3 and the equilibrium field- 
induced state is reached at tt = 30 x 10 4 . The correspond- 
ing free-energy is shown in (a), top panel. Then the field 
is switched off and the system relaxes to a new metastable 
state (tg = 80 x 10 4 ), ((b) in the bottom panel) character- 
ized by a new free energy, slightly higher than the bulk BPI 
value. An electric field along y is then switched on and off at 
ii2 = 120 x 10 4 ((c) and (d) in the top panel). Strikingly, the 
relaxation process brings the system back to the equilibrium 
state without field. 



ity, and the switching dynamics of the defect network is 
reported in Fig. [3] When the electric field is switched on 
along the x direction {t\ — £4), the defects start to deform, 
bend and form a spectacular double-helix pattern in the 
bulk, accompanied by columnar defects at the wall. This 
structure is induced and stabilised by the electric field, 
yet it only suffers minor rearrangements as the field is 
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switched off. This double-helical state has already been 
observed in Ref. |15) (a) in the absence of a field - it could 
be selected for instance by tuning the device thickness. 
When we switch on the field again, now along the y direc- 
tion (see snapshots tg — tie), the bent double-helix moves 
towards the walls (tg — t w ) where it joins up with the 
the columnar disclinations, to form a transient intercon- 
nected structure (in) which later on relaxes to an array 
of twisted ring defects in the bulk, accompanied again 
by straigth disclinations, this time at right angles with 
the one associated with the double helical state (£12). 
Once more the field-induced rings do not change much 
upon switching off (t\2 — iig)- We note that a state sim- 
ilar to these ring defects has recently been observed in 
Ref. |15)(b). however in a system with homogeneous an- 
choring at the walls. Our simulations therefore show that 
the number of possible defect networks in these thin BP 
cells are even more than previously thought. In order to 
behave as a switchable bistable device, our system still 
needs to be able to reconstruct the double-helical state 
starting from the ring defect structure. This is indeed 
what happens if the field is switched on and off along the 
x direction (as in the first cycle). The ring defects ini- 
tially enlarge and touch the walls where disclination lines 
reform (£17 — iis). Thereafter the double-helix pattern re- 
constructs in the bulk (tig) and later on stabilises. Thus 
our thin BP sample with homeotropic anchoring behaves 
as a bistable device, and the two switchable states are 
the double helical and twisted ring structures - unlike 
the surface memory device none of these equals the zero- 
field stable state. 

In summary, we have investigated the hydrodynamics 
of a cubic BPI cell and proposed two systems, differing in 
the boundary conditions, which can switch between two 
bistable disclination structures under the action of an ap- 
plied field. In the first one we have fixed boundary condi- 
tions at both walls, while in the second one we considered 
homeotropic anchoring. In the FBCs cell, the two switch- 
able states are the zero-field BPI and another metastable 
defect network. The device with homeotropic anchoring, 
on the other hand, switches between two states which are 
metastable in the absence of a field. In one of these states 
the defects form a double helix, whereas in the other one 
they arrange as an array of twisted rings - these config- 
urations are alternatively observed by varying the direc- 
tion and the magnitude of the electric field. The bistable 
behaviour of this device is promising for technological 
applications to energy-saving devices because, in prin- 
ciple, the required anchoring could easily be achieved. 
The physical reason behind the viability of BP bistable 
devices may be the glassy free energy profile which ad- 
mits a number of competing metastable minima, each of 
which is potentially suited for their design, provided a 
suitable kinetic route to achieve reversible switching is 
discovered. In this respect, one may expect that even 
more BP bistable systems may be found in the future. 
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FIG. 3: Defect evolution with homeotropic anchoring un- 
der an applied electric field. Parameters are the same as in 
Fig. [2] The electric field (~ 4.5 V/fim, i.e. with e 2 ~ 0.06) 
is switched on along x at time t = 10 X 10 4 (Fig. [TJd for the 
steady state) and the evolution dynamics is reported in the 
first row. A steady state double-helix defect is attained at 
time ti = 98.7 x 10 4 . When the field is switched off the 
defect structure is almost unaltered (steady state reached 
at time t 8 = 117.7 x 10 4 ). Then field (~ 4 V/fim, i.e. 
with e 2 ~ 0.05) is then switched on along y. The double- 
helix disclination transforms into a steady state with twisted 
ring defects (ti2 = 155.6 x 10 4 ) stable after switching off 
(tie = 194.2 x 10 4 ). Bistability requires a route back to the 
double helix state: this is provided by a further application 
of a field along x (steady state at £20 = 223.4 x 10 4 ). 
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